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Abstract 

We consider (AlAs) n /(GaAs) n superlattices with random thickness fluctua¬ 
tions An around the nominal period n. Using three-dimensional pseudopoten¬ 
tial plane-wave band theory, we show that (i) any amount A n/n of thickness 
fluctuations leads to band-edge wavefunction localization, (ii) for small A n/n 
the SL band gap is pinned at the gap level produced by a single layer with 
“wrong” thickness n + An, (iii) the bound states due to monolayer thickness 
fluctuations lead to significant band-gap reductions, e.g., in n = 2,4,6, and 
10 monolayer SL’s the reductions are 166, 67, 29, and 14 meV for (111) SL’s, 
and 133, 64, 36, and 27 meV for (001) SL’s, (iv) (001) AlAs/GaAs SL’s with 
monolayer thickness fluctuations have a direct band gap, while the ideal (001) 

SL’s are indirect for n < 4. 
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The electronic structure and quantum-confinement effects in semiconductor superlattices 
are usually modeled by assuming an ideal structure, i.e., that the interfaces are atomically 
abrupt and that the individual layer thicknesses remain constant throughout the superlat¬ 
tice [[L],f2|. Actual heterostructures, however, often deviate from ideality in two ways: (i) 
lateral imperfections in the (x,y) plane, such as intermixed ||], stepped [f|, or islanded 0 
interfaces, and (ii) vertical ( z direction), discrete thickness fluctuations around its nominal 
value, even though the interfaces could remain reasonably flat and atomically abrupt in 
the ( x,y ) plane ||. It is likely that both types of interfacial imperfections coexist in many 
samples []7|]. In case (i) the translational symmetry is broken in the substrate plane (x,y), 
and the concentration profile along the growth direction is continuous (“graded” or “in¬ 
termixed” interfaces), while in case (ii) a discrete, rectangular-shaped concentration profile 
exists, but the superlattice translational symmetry is broken along the growth direction z. 
Case (i) has been modeled theoretically by assuming graded (rather than square) potential 
wells |8|], or by considering supercells with atomic swaps across the interfaces [|9|. In this 
paper we consider layer-thickness fluctuations [case (ii)], i.e., ( A) n /(G) n superlattices (SL) 
with nominal layer thicknesses of n monolayers (ML) of material A = AlAs and n ML of 
material G = GaAs, but where each layer is occasionally thinner or thicker than its intended 
thickness n. We use a three-dimensional pseudopotential band-structure description within 
a highly flexible plane-wave basis [|K|, rather than one-dimensional effective-mass models 
m , or one-dimensional |12| or three-dimensional [|L|] tight-binding models. 

Superlattices with layer-thickness fluctuations are described by supercells containing sev¬ 
eral hundreds of atoms |T4|. The novel empirical pseudopotentials used here [|T0|] have been 
tested extensively for AlAs/GaAs bulk materials, short-period superlattices, and random al¬ 
loys. The results [[U]] compare well with experiment and with state-of-the-art, self-consistent 
pseudopotential calculations, however, without suffering from the band-gap underestimation 


of the local-density approximation [TIJ. For periods n < 20 we consider single monolayer 
fluctuations , so the layer thicknesses are in the set {n — 1 ,n,n + 1}, while for n > 20 we 
consider a fixed fraction A n/n of layer-thickness fluctuations An. 
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We find that: (i) any amount A n/n of thickness fluctuations leads to band-edge wave- 
function localization, (ii) for small A n/n the SL band gap is pinned at the gap level produced 
by a single layer with “wrong” thickness v! 7 ^ n (a “chain mutation”), (iii) the bound states 
due to monolayer thickness fluctuations lead to band-gap reductions that monotonically 
decrease with increasing n. These fluctuation-induced bound states will photoluminesce at 
energies below the “intrinsic” absorption edge, (iv) (001) AlAs/GaAs SL’s with monolayer 
thickness fluctuations have a direct band gap, while the ideal (001) SL’s are indirect for 
n < 4. 

Consider first, by way of reference, the band structure of ideal AlAs/GaAs SL’s with 
layers oriented along (111) or (001) (Fig. |l]). In agreement with previous theoretical studies 
0, we find that (i) the ideal (111) SL has a direct band gap for all n values, since the 
conduction-band minimum (CBM) is the T-folded r lc (r lc ) state. This SL has a “type-I” 
band arrangement with both the highest valence and the lowest conduction state localized 
on the GaAs layers, (ii) The second conduction band at T is folded from the zincblende T-L 2 
bands; for small n the pseudodirect f lc (L lc ) state mixes strongly with the direct T lc (T lc ) 
state. The mixing, and thus the level repulsion, shows odd-even oscillations for small n 
(reflecting localization of repelling states on the same or on either sublattice [^J). Note that 
the one-dimensional effective-mass model (dashed lines [l5j in Fig. |l|) completely misses the 
strong non-monotonic variations of SL energy levels with layer thickness. 

The situation is very different for (001) oriented ideal SL’s. The prominent properties 
apparent in Fig. |l](b) are: (i) the n = 1 SL has an indirect band gap at the L-folded point R 
[l].[2j]; (h) for n < 4, the lateral X xy valleys (folded to M) and the X z valley (folded to F) are 
nearly degenerate [16|]; (iii) for 1 < n < 8 the pseudodirect, AlAs-like fi c (AL) state is below 
the direct, GaAs-like r lc (r lc ) state, thus the SL is type II; for n > 8 , however, r lc (r lc ) 
is lower, so the system is type I (experimentally, the type-II/type-I crossover is found at 

n ~ 11 0 ); 

We now allow the layer thicknesses in the n X n SL’s to fluctuate around the ideal value 
n by ±1 ML. The growth sequence is now defined in terms of a distribution function p(n'), 
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which we assume to be uncorrelated and symmetric around the nominal thickness n. We 
define the relative frequency R of monolayer fluctuations by 

p(n ± 1) = Rp(n). (1) 

Because the distribution p(n') is normalized, we can write p(n) = 1/(1 + 2 R) and p(n± 1) = 
R/( 1 + 2 R). For the ideal n x n superlattice R — 0, and p{n') = 8{n' — n), whereas for 
R — 1 the layer thicknesses {n — 1 ,n,n + 1} occur with equal probability p = \. A single 
chain mutation in a finite superlattice of length N monolayers corresponds to R ~ which 
will denoted here as the R —> 0 limit (to be distinguished from R = 0 with no mutations). 
To simulate the lack of periodicity along the growth direction, we have used supercells of 
a total length N between 100 and 1000 ML, and repeated these supercells periodically. A 
particular growth sequence was created by a random number generator, specifying n, N and 
R. 

The band-edge energies of (AlAs) n /(GaAs) n SL’s with one-monolayer thickness fluctua¬ 
tions about n are plotted in Fig. |2| relative to the band edges of the ideal SL’s (the energy 
zero). We see that: (i) thickness fluctuations create both hole (Aeh) and electron (Ae e ) 
bound states for any degree of thickness fluctuation (0 < R < 1), (ii) the band-gap re¬ 
ductions AE S = A^h + Ae e decay with n, and have a definite dependence on superlattice 
direction; they are 166, 67, 29, and 14 meV for n =2,4,6, and 10 in the (111) direction, and 
133, 64, 36, and 27 meV in the (001) direction, respectively, (iii) the dilute limit of a single 
chain mutation already produces a finite gap reduction A E g (R —» 0), (iv) A E g (R —> 0) 
merges with A E g (R = 1) at n > 6, at which point the gap reduction becomes independent 
(“band-gap pinning”) of the number of chain mutations. 

The appearance of gap levels in SL’s with one-monolayer thickness fluctuations is accom¬ 
panied by wavefunction localization. For example, inspection of the CBM wavefunction of 
an n — 6 (111) superlattice with random ±1 monolayer fluctuations [Fig. |^(a)] reveals that 
it is localized on ~4 GaAs wells, with minimal amplitude in the AlAs barriers and maximal 
amplitude on the two neighboring mutated (7-ML) GaAs wells (“twin” fluctuation denoted 
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by bold arrows). The CBM thus resembles a bound state in a coupled double quantum 
well. The hole wavefunction at the valence-band maximum (VBM) is likewise localized on 
a number of mutated, 7-ML GaAs wells [Fig. ^j(a)]; in contrast to the CBM, however, the 
reason for the multi-well pattern of the VBM wavefunction is that these states are in fact 
decoupled, quantum-well confined states, which are degenerate in energy within the accu¬ 
racy of our calculation (<0.1 meV). A typical hole and electron wavefunction localized on 
an isolated (GaAs )7 mutation in an otherwise ideal 6 x 6 (111) SL are shown in Fig. ||(b). 
We see that the hole wavefunction of an isolated mutation (R —> 0) resembles that of the 
concentrated (R = 1) mutations [Fig. |3|(a)], and its binding energy Ae^R —► 0) = 11 meV 
equals the value at R = 1. At the CBM, the larger penetration of the wavefunction into 
neighboring GaAs wells can produce deeper gap states, and consequently pinning occurs at 
a larger n ( n p ~ 10 ) than for hole states. 

Experimentally, the fluctuation-induced localized bound states should be observable as 
photoluminescence centers whose energy is below the absorption edge of the underlying 
“ideal” SL structure. This photoluminescence will lack phonon lines, because the optical 
transitions are direct in the planar Brillouin zone (the transverse wave vector k± is still a 
good quantum number), and because the k z selection rule is relaxed by vertical disorder. 
Indeed, while calculations f 2 | on ideal (111) (AlAs) n /(GaAs) ri SL’s predicted a direct band 
gap with a type-I band arrangement, Cingolani et al. 13 noted a ~100 meV red shift of 
the photoluminescence at 1.80 eV relative to the absorption in (AlAs) 6 /(GaAs ) 6 (111) SL’s, 
interpreting this as reflecting a type II band arrangement. However, since they noted that 
their SL had a ±1 ML period uncertainty, it is possible that the red shifted photolumines¬ 
cence originates from thickness-fluctuation bound states. Our calculated band gap of the 
n — 6 superlattice with ±1 ML thickness fluctuations is 1.78 eV for R = 1, and 1.80 eV for 
R —> 0, close to their observed photoluminescence peak position (1.80 eV) |17| , [T8 . 


Figure |2j shows that the bound states of isolated mutations (R —> 0) merge with those 
of concentrated layer-thickness fluctuations (R —> 1) at some “pinning period” n p . In what 
follows we discuss the (i) n < n p and (ii) n > n p regimes. I 11 the short-period regime 
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(n < rip), the band-gap reduction depends on R. This regime coincides with the region in 
Fig. 0 where the band gaps of the ideal SL’s have a complex n-dependence, showing strong 
non-effective-mass behavior. The (001) (AlAs) 2 /(GaAs )2 SL with monolayer fluctuations 
is in fact identical to the intentionally disordered SL grown by Sasaki et al. In that 

structure, A 2 and G 2 layers are randomly replaced by A^,A 3 ,Gi, and G 3 layers. We find 
the following changes in the band structure when the layer thicknesses fluctuate by ±1 ML: 
(i) For the n = 2 (001) SL we obtain A s e (R = 1) = 22, 81, and 171 meV at M, f, and X in 
the planar Brillouin zone. Since the level shift Ae e at T exceeds the one at M by ~60 meV, 
layer-thickness fluctuations transform the indirect 2x2, ideal superlattice into a direct-gap 


material p0[. The large binding energy at X originates from the increased level repulsion 
of the folded L lc states when the translational and rotational symmetry of the ideal n — 2 
superlattice is broken. This level repulsion is larger for odd n than for even n, and it leads 
to an L-like, indirect CBM for n — 1 [Fig. [|(b)]. (ii) The n = 2 (001) SL is direct even in 
the R —> 0 limit of isolated chain mutations, (iii) Numerous electron and hole states are 
localized on the same spatial region along z, hence the (001) SL’s with thickness fluctuations 
will exhibit type-I, rather than type-II characteristics as ideal n < 10 SL’s do. (iv) While for 
n > n p only mutated, wider GaAs wells bind a carrier (see below), in the short-period (001) 
case even an (AlAs) n+ i mutation binds an electron. In fact, the AlAs-like bound electron 
lies deeper in the gap than the GaAs-like bound electron [see the two dotted lines near the 
CBM in Fig. @(b)]. 

We next discuss the SL properties in the regime of band-gap pinning n > n p . At this 
point the band-gap reduction is pinned at the value 


A E g (R) = lim AE g (R) = Ae e + Ae h , (2) 

where Ae e (Af^) is the electron (hole) binding energy of an isolated ( R —> 0) layer mutation. 
Qualitatively, Eq. (^) can be understood in terms of a one-dimensional effective-mass picture, 
where the SL is modeled by a Kronig-Penney model, with the usual boundary conditions of 
a continuous envelope function F(z) and current A^F\z), where the effective mass m* is 
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different in the well and barrier material. Each of the (n + 1)-ML mutations gives rise to a 
bound state below the band edge of the nxnSL BED- The gray lines in Fig. ^indicate that 
this simple picture is expected to agree quantitatively with our pseudopotential calculation 
for n ^ 10 (note, however, the large discrepancy with our three-dimensional calculation at 
smaller n’s). For very large n, when the quantum wells are completely decoupled, the SL 
energy spectrum is simply that of degenerate single quantum wells. Hence the extra binding 
energy of an (n + An)-ML mutation approaches asymptotically 

^/\79 

Ae c = £ 0 (n) -£ 0 (n + An) «- e 0 (n), (3) 

n 

where £o(n ) is the ground-state energy of a carrier with mass rn* in a n-ML wide quantum 
well, which scales like —^ for large n. Using a fixed A n/n = 10% we obtain from the 
first equality of Eq. (||) Ae e = 10.0,2.4 and 0.7 meV for n = 20,50 and 100 in the (111) 
SL [the last equality of Eq. (^|) gives 14.3, 3.0 and 0.8 meV, respectively]. The band-gap 
reduction for a given A n/n is obtained by inserting A^h and Ae e from Eq. (||) in Eq. (^|). 
The degeneracy of the gap level Ae c is equal to the number M n+ ^ n of (n + An)-ML well 
mutations, which, in the case of An = 1 [see Eq. ([!])], is given by 


M, 


R N 


n+1 


1 + 2 R 2 n 


( 4 ) 


Using Eqs. (|3|) and (f|) one can predict the band-gap reduction and number of localized 
states in the multiple quantum-well regime. 

In summary, we have shown that discrete layer-thickness fluctuations in (AlAs) n /(GaAs) n 
superlattices lead to (i) localization of the wavefunctions near the band edges, (ii) a reduced 
band gap, which is pinned at the value corresponding to an isolated layer-thickness fluctu¬ 
ation for n > 7i p , and (iii) a crossover to a direct band gap in the case of short-period (001) 
SL’s. 
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FIGURES 


FIG. 1. Energy levels of ideal (AlAs) n /(GaAs) n superlattices along (a) (111) and (b) (001), 
as a function of period n. The bulk levels in the middle column are reached asymptotically as 
n —> oo. ri 5 V (AlAs) is 500 meV below the GaAs VBM. The SL states are denoted with an overbar, 
while their parent zincblende states T, X , and L are given parentheses. The dashed lines are 
obtained from a one-band envelope-function model using the effective masses and band offset from 


our pseudopotential calculation 115]. 


FIG. 2. Gap levels in the presence of one-monolayer thickness fluctuations in (AlAs) n /(GaAs) n 
superlattices along (a) (111) and (b) (001), as a function of period n. Energies are measured with 
respect to the band extrema of the ideal nxn SL (see Fig. |l|). One-band envelope-function results 
for a (n + 1)-ML mutation embedded in the nxn SL are indicated by a dashed gray line. R = 1 
and R —> 0 denote, respectively, the concentrated and dilute limit of chain mutations [Eq. #]• 


FIG. 3. Planar averages of wavefunctions squared of the CBM and VBM in the 
(AlAs)e/(GaAs)6 SL along (111) with ±1 layer-thickness fluctuations. Hole wavefunctions are 
plotted in the negative direction, with a small offset for clarity, (a) Concentrated limit (R = 1), 
(b) dilute limit (R —► 0), i.e., a single (GaAs )7 mutation embedded in a 6 x 6 SL host. The rectan¬ 
gular lines show the growth sequence of the SL, with GaAs layers represented by wells, and AlAs 
layers represented by barriers, respectively. The vertical arrows in (a) indicate the 7 ML thick, 
“mutated” wells. 
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